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We show that the recently proposed doped carrier Hamiltonian formulation of the t-J model 
should be complemented with the constraint that projects out the unphysical states. With this new 
important ingredient, the previously used and seemingly different spin-fermion representations of 
the t-J model are shown to be gauge related to each other. This new constraint can be treated 
in a controlled way close to half-filling suggesting that the doped carrier representation provides 
an appropriate theoretical framework to address the t-J model in this region. This constraint also 
suggests that the t-J model can be mapped onto a Kondo-Heisenberg lattice model. Such a mapping 
highlights important physical similarities between the quasi two-dimensional heavy fermions and the 
high-Tc superconductors. Finally we discuss the physical implications of our model representation 
relating in particular the small versus large Fermi surface crossover to the closure of the lattice spin 
gap. 
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I. INTRODUCTION 

The high-Tc superconductors continue to be a puzzle 
to most researchers in the field. Notably, in the un- 
derdoped regime, the cuprates display highly anomalous 
physical properties. Above the superconducting temper- 
ature, Tc, those lightly doped materials are characterized 
by a spin gap and, by what seemed until recently, discon- 
nected Fermi arcs around preferential directions in mo- 
menta spaceii^. This phase is widely referred to as the 
pseudogap state. That such disconnected arcs are indeed 
integrated into one coherent Fermi surface was recently 
demonstrated by the experiment of Doiron-Leyraud et 
al2, unfolding the quantum oscillations in the electrical 
resistance of YBa2Cu306.5. The pseudogap phase is in 
this way associated with a small Fermi surface as opposed 
to the overdoped cuprates which exhibit large Fermi sur- 
faces instead^. The evolution of the cuprate Fermi sur- 
face as a function of doping has been monitored by sev- 
eral ARPES experiments^'*'. Such a process has also been 
the object of investigation in a number of renormalization 
group calculations^. Despite that, it is fair to say that 
a more qualitative physical understanding of this Fermi 
surface crossover in the cuprates is still lacking. 

Ever since Anderson's suggestion that the high-Tc 
cuprates are driven by strong electron correlations^, the 
Hubbard model, or its large U version, the t-J model, 
have been the center of attention of theorists. Re- 
cently, two new representations in terms of dopant par- 
ticles have been proposed^'^'^'^^ for the t-J model. In 
the Hamiltonian approach put forward by Ribeiro and 
Wen^iiS, the original projected electron operators are re- 
placed by spin-1/2 objects (called "lattice spins") while 



the dopant particles are represented by fermions with 
spin-1/2 (called "dopons"). In the path-integral repre- 
sentation, the spin-dopon partition function of the t-J 
model has been used to formulate a resonating valence 
bond mean-field theory to describe the superconducting 
phase in the cobaltatesii. The motivation behind these 
approaches is clear: since the concentration of dopons 
is small close to half- filling, the no double occupancy 
(NDO) constraint for the dopons can be safely relaxed 
in that regime. 

However, in the description of electrons in terms of 
spins and dopons, special care should be taken to avoid 
the inclusion of unphysical states. As pointed out by 
Ferraz, Kochetov and Uchoai^,, a constraint should be 
imposed to eliminate the unphysical states. Within the 
doped carrier representation, the conventional no double 
occupancy constraint for the lattice electrons reasserts 
itself as the constraint to eliminate the unphysical states 
in the enlarged spin-dopon Hilbert space. Basically, it 
reflects the completeness relation of the physical Hilbert 
space for the t-J model. 

It is important to note that in the Hamiltonian doped 
carrier representation of the t-J model the NDO con- 
straint takes on a form which is in a sense dual to that 
of the original one in the physical lattice electron rep- 
resentation. The original local NDO constraint for the 
physical lattice electrons, '^^cla-Cia < 1, cannot be re- 
laxed at half- filling (in fact, right at half-filling it reads 

J^a'^la'^i'^ — 1)' but it can be dropped safely at a large 
enough doping. In contrast, the new spin-dopon con- 
straint can be treated in a controlled way close to half- 
filling (right at half-filling this constraint becomes a triv- 
ial identity), whereas it cannot be safely relaxed in the 



2 



overdoped region. Since the underdoped region of the 
cuprate phase diagram is of primary interest as was 
stressed earher on, this new doped carrier formulation 
of the t- J model, accompanied with the spin-dopon con- 
straint, sounds quite interesting and appealing. 

Within the spin-dopon representation, the strong elec- 
tron correlation manifests itself in the constraint to ex- 
clude the unphysical states. It turns out that this con- 
straint can be treated reliably at the mean-field level close 
to half-filing. This occurs at the expense of having more 
complicated form of the t-J Hamiltonian in this new rep- 
resentation. 

The present paper demonstrates the significance of 
the constraint that excludes the unphysical states in the 
doped carrier approach from both the physical as well 
as the computational points of view. Specifically, we 
show that there is a certain gauge ambiguity in choos- 
ing a particular form of the unprojected t-J dopon-spin 
Hamiltonian in the enlarged Hilbert space. This ambigu- 
ity is related to a redundancy of the representation of the 
projected electron operators in terms of the lattice spins 
and dopons. This ambiguity is removed by projecting 
the gauge-dependent Hamiltonian onto the physical sub- 
space: upon projection all the gauge-related Hamiltoni- 
ans result in one and the same physical representation. In 
this way we show that the previously used and seemingly 
different spin-fermion representations of the t-J model 
are in fact related to each other by this projection. 

The spin-dopon constraint corresponds to a Kondo- 
type interaction between the lattice spins and the dopons. 
The dopons play the role of the "conduction electrons" 
while the lattice spins play the role of the "localized mag- 
netic moments" . This approach allows for a mapping of 
the original t-J model onto a Kondo-Heisenberg lattice 
model which indicates a strong relationship between the 
physics of high-Tc materials and that of some of the quasi 
two-dimensional heavy fermion systems (see sections III- 
IV). We also show that this constraint, in contrast to the 
standard NDO condition, can be treated in a controlled 
way within mean field approximation. 

Using our approach, the pseudogap phase can be sim- 
ply interpreted as a state in which the lattice spins are 
paired, and the low density dopons are the quasiparti- 
cles solely responsible to generate the small Fermi sur- 
face. If the dopon concentration increases so does the 
lattice spin-dopon coupling. As a result, in the normal 
state above Tc there is a critical dopon concentration 
at which this coupling starts to dominate over the cor- 
responding lattice spins self-interaction finally breaking 
the remaining spin pairs. When this takes place, both 
spins and dopons integrate the large Fermi surface and 
the physical system crosses over to a metallic phase with 
no pseudogap behavior. 



II. DOPED CARRIER REPRESENTATION 

The t-J model hamiltonian in terms of the Gutzwiller 
projected lattice electron operators takes the form 

-fftJ = - ^ Uj4a^ja + • Qj - ^n,nj), (1) 

ija ij 

where Cia = Ci(j(l — ni_^a) is the projected electron op- 
erator (to avoid the on-site double occuppancy), Qi = 
Sct (T ^L''CT(TCiCT, = 3/4, is the electron spin operator 
and hi = -|- n^j — 2ni|ni|. 

Following Ribeiro and Wen, we consider an enlarged 
Hilbert space for each site i. This enlarged space is char- 
acterized by the state vectors |cra) with a =t, i labelling 
the spin projection and a = 0, t,i labelling the dopon 
state (double occupancy is not allowed). In this way the 
enlarged Hilbert space becomes 

W"^"' = {| T o)„ I i o)„ I Ti)„ I it)., I TT)., I ii) J (2) 

while in the original Hilbert space we can either have one 
electron with spin cr or a vacancy: 



w = {|T).,U).,|o),} 



(3) 



The following mapping between the two spaces is then 
defined: 



T)^^ I TO). 



li). U 0). 



ITi).-UT). 

V2 



(4a) 



(4b) 



(4c) 



The remaining states in the enlarged Hilbert space, 
ITi).4^J-T). ^ I 11^.^ I 11^. g^j.g unphysical and should be 

excluded in practical calculations. In this mapping, a 
vacancy in the electronic system corresponds to a singlet 
pair of a lattice spin and a dopon whereas the presence 
of an electron is related to the absence of a dopon. 

Let us call Si the operator associated with the lattice 
spin on site i and dj^ the creation operator for dopons. 
Then under this mapping it is not difficult to find rela- 
tions between the original projected electronic operators 
and the new operators, such as: 



Ci] 



1 

72 



(5) 



where diu = dia{l — d\ _^di,^^) is a projected dopon op- 
erator. Although Si and di^^ act in the whole enlarged 
Hilbert space, specific combinations like the one given in 
Eq. (O nullifies the unphysical states. It is then possible 
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to write the original t- J Hamiltonian in terms of the new 
operators in such way that it vanishes when acting on 
the unphysical states — . The reason for this construc- 
tion is now obvious: close to half-filling there is a small 
amount of dopons in the system and as a result the local 
constraint of no double dopon occupancy can be safely 
dropped. 

However, as soon as some mean-field approximations 
are performed, the unphysical states reappear in the the- 
ory regardless of whether projected or unprojected dopon 
operators are used. A constraint should therefore be im- 
posed in order to eliminate the unphysical states. This 
constraint was proposed in and it reads 



-I -> 3 ~ ~ 
Si ■ Ml + -(d||diT 







(6) 



where Ali = ^ dl^Ta-a-dia- is the spin associated with 

the dopon. Acting on the physical states, gives zero 
while T i\unphys)i = \unphys)i. 

Since for 1/2 spin 5f = 3/4 and = (3/4)nf, this 
constraint can also be recast into the form 

Jt - 3/4(1 - nf) = 0. 

Here Ji — Mi + Si is the total spin on each lattice site. 
In this form the constraint tells us that the on-site total 
spin can be cither j — (vacancy) or j ~ 1/2 (lattice 
spin) . 

There is a yet another representation of this constraint, 
namely 



(7) 



where d\ — — /ji'^ji creates an on-site vacancy. 

Here we represent the spin degrees of freedom in terms 
of chargeless fermions (spinons), 5^ — f-ffi, fi = 

ifli, fliY, fifi = 1- In this form the constraint Y.i'^ = 
tells us that the total number of vacancies must be equal 
to the total number of dopons^^. 

Right at half-filling {nf = 0), the spin-dopon con- 
straint becomes a trivial identity, = 0. It can thus be 
safely treated close to half-filling at a mean-field level. In 
contrast, in the overdopcd regime {iif w 1) it reduces to 
the equation, 5*^ - Mi w —3/4. Exactly the opposite situa- 
tion is realized for the lattice electron NDO constraint: it 
cannot be relaxed close to half-filling but it can be totally 
ignored for a low enough density of the lattice electrons. 
It is precisely because of this that the new dopon-spin t-J 
model representation becomes indeed quite appealing. 

It is convenient to define the operator Pi — 1 — T^. 
Since Tf = Ti it is clear that Vi is a projection operator 
which eliminates unphysical states on the site ri^. We 
show now that by means of this projection operator we 
can indeed establish alternative ways of expressing the 
electron operators in terms of the lattice spins and do- 
pons. Let us consider for example the action of dii on the 
physical states - it is simply d^J f 0)i = 0, dij i 0)i = 



and Jixll^ 



V2 



Clearly, in the physical sub- 



space this is equivalent to the action of cJ|/\/2. Thus we 
can readily write: 



(8) 



Alternatively we can use the projection operator Vi to 
perform explicit calculations, resulting in 



V2-Pid,iV, 
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S+d.. 



•it 



which combined with the adjoint of Eq. ^ results again 
in our Eq. ([8]). In the same way one can show that 



^V2Vrd,^Vi. 



(9) 



Consider now the action of S^di-^ on the physical 
states, namely, S^di^\ t 0)i = 0, S^di^^ \ J, 0)i = and 
^r+J^^lUh^il). ^ -i^. Comparing with the action of 

c|^|/\/2, or carrying out the explicit calculations with the 
V projection we get yet another operator identity: 



(10) 



whereas from the action of S^ d^ we immediately find 



(11) 



Note that a few attempts have also been made at de- 
coupling the physical electron as the spinful fermion and 
spinon. However, they can only be justified within some 
approximation scheme (see, e.g., ref.— ). In contrast, here 
we display the exact form of the spin-fermion decoupling 
of the projected lattice electron operator in Eqs. (10)- 
(11). In the next section we show that the seemingly 
different representations for one and the same lattice elec- 
tron operator (e.g., given by Eqs. (9) and (11)) are in fact 
identical since the V projections of gauge related objects 
coincide with each other. 



III. LOCAL GAUGE SYMMETRY 

Let us now define the global projection operator V ~ 
Yli'Pi- The t-J Hamiltonian can be generally written 
as Htj = VHV, where H acts on the enlarged Hilbert 
space. From the discussion above it is clear that there are 
different choices of H that after being projected onto the 
physical subspace result in the t-J model Hamiltonian. 
This is related to the gauge ambiguity of the unprojected 
H with respect to the local U(l) gauge transformations 
generated by the constraint Ti. Under this transforma- 
tion. 



i) = exp{-i6iTi)di^iexp{i9,T,) 



4 



S-di. + {- + S!)d^, 



(e^"' - 1) + 



(12) 



S+d, 



ii 



(e'^--l) + dx,. (13) 



The explicit form of the gauge dependent operator Si{9) 
is given by 



Sihi - 



+ 



cos6li - 1) + Si 



(14) 



Here d„., = d„i[6 = 0), Si = S^{e = 0). The U(l) gauge 
symmetry is reahzed on the spin-dopon multiplet in a 
nontrivial way. Note also that the total on-site electron 
spin operator Ji — Si + Mi as well as the dopon number 
operator are gauge invariant quantities, Ji{9) = Ji{6 = 
0), nm = ni{9 = 0). 

The ■p-projection of the different gauge equivalent op- 
erators results in the same gauge invariant representa- 
tion. For instance, for any polynomial on-site operator 
fi = Mdai,Si) we get 

-PifiiO)A = exp(-z0,t,)/. exp(z0,t,)Pi 



= r, [i + r^ie'"' - 1) J /, [ 1 + v,ie-''' - i)J = r,m. 

In this way the representation (11) follows from the equa- 
tions 

Adt^A = 'P^d1i,{e = Tr)?^, = -Vr{S^dli)Vr 



= 2 - S^ dii 

Notice that in general 



'PY.*.,dlmd,Ae,)'P 

ija 



ija 



It can also be checked by an explicit computation that the 
projected electron operators, e.g., the ones given by Eqs. 
([51 ITT|) are all gauge invariant objects. In analogy with 
the gauge theories from quantum field theory, we can say 
that essentially the choice of a given representation for 
H corresponds to the fixing of a particular gauge. 

The existence of the local gauge symmetry reflects 
a degree of redundancy in the parametrization of the 
Gutzwiller projected lattice electrons in terms of the lat- 
tice spins and dopons, as displayed in our Eq. ([5]). In 
principle, one can formulate a mean-field theory in this 
representation that respects that local U(l) gauge sym- 
metry. In this way one arrives at a local gauge theory 
that describes quantum fluctuations around the mean- 
field solution. 



It should be kept in mind that the lattice spins and 
dopons are, in general, not gauge invariant and couple to 
the gauge field. Because of this the dopons and lattice 
spins, away from half-filling, do not represent real excita- 
tions and they are introduced as an intermediate step to 
calculate the physical (gauge-invariant) quantities such 
as given, e.g., by Eqs. (O, (8)-(ll). Note however that, 
right at half- filling. Si ^ Si. In other words the lat- 
tice spins represent, in this limit, real excitations. It is 
therefore natural to assume that close to half-filling the 
lattice spins and dopons can be viewed as well-defined 
excitations weakly coupled to the gauge field. This indi- 
cates that the mean-field spin-dopon theorj>ii is presum- 
ably stable close to half-filling with respect to quantum 
gauge fluctuations. However, an explicit estimation of 
the strength of the gauge interaction can be made only 
after the full gauge theory is derived. The explicit form of 
that theory is not still available because of a rather com- 
plicated form of the U(l) group action on the spin-dopon 
multipletes (see Eqs. (jl2m4p V 

To make this point more clear, let us contrast the prop- 
erties of the gauge symmetry generated by Ti with those 
of the U(l) local gauge symmetry generated by the stan- 
dard NDO constraint in the frequently used slave-boson 
representation for the lattice electron operators. With 
this formalism the projected electron operator is written 
as 



with the NDO condition 



t 



fib. 



(15) 



where fia is a fermion (spinon) operator and hi is a 
slave-boson (holon) operator. Conservation of the gauge 
charge Qi can be derived by the Noether theorem start- 
ing from the local U(l) gauge transformatiouji^ 

.Ua ^ e''*/.., h - e^'-h, (16) 



which leaves the physical electron operator (|T5)) intact. 
This U(l) local gauge symmetry takes care of the redun- 
dancy of the parametrization (|15p . In contrast to the 
spin-dopon charge T^, the operator Qi does not vanish 
at half-fiUing, Qi = Y.afLfi'^- This indicates 

that the spinous are strongly coupled to the bare gauge 
field close to half filling. However, the auxiliary gauge 
field can in general acquire nontrivial dynamics at low 
energies, effectively moving the model to a weak cou- 
pling regime. Therefore, to judge whether confinement 
or deconfinement of slave particles really occurs in the 
physical low-energy excitations, one must explicitly in- 
vestigate the gauge dynamics in the low-energy regime^. 

Note also that in contrast with the NDO constraint for 
the lattice electrons, the set of local constraints Ti = 
(one for each lattice site) can be replaced by the global 
condition T = J^i '^i ~ 0- The reason for this simplifica- 
tion is that the unphysical states appear as the degener- 
ate eigenvectors of T.^ with an eigenvalue, 1. Therefore, 
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acting on an unphysical state, T simply produces the 
same state multiplied by a positive number. Acting on 
a physical state, T always gives zero. Note, however, 
that Ti involves a quartic power of interacting dopons 
and spins. The standard NDO constraint appears as a 
quadratic form of the electron operators. We will com- 
ment on this point further at the end of the paper. 

Summing up all the above, we can write down the exact 
form of the t-J Hamiltonian in the spin-dopon represen- 
tation as 



-J^^((5,+Al.)(5,+Al,)--(l-n^)(l-n^^))7^ (17) 



Within the path-integral approach this representation 
has been used inii to obtain the mean-field Tc phase 
diagram for the cobaltates. The exact spin-dopon path- 
integral representation of the t-J partition function given 
irtii is written down in terms of the classical fermion am- 
plitudes which are related to the dopon amplitudes in 
the following way, V« = -^i^-2SfT)df = V2{VAV^Y'- 
This relation holds provided Tf ~ 0. Note that in this 
case (nfy' = {fifyK Within that path-integral approach 
the constraint to exclude the unphysical states reasserts 
itself in the form of the SU(2) invariant site product of 
the delta-functions that singles out the physical subspace. 

Since Q,^r,{S^ + M,)r^ = 5,(l-d|^d,<,), and 
d\^di„)Vi — — dl^dia), we can rewrite Eq.([T7|) in the 
form 

Htj=vY.'^U,dld,,V 

ija 



(18) 



Note the important factor of 2 in front of the t-term in 
these formulas. It originates from the fact that the va- 
cancies are represented in this theory by the spin-dopon 
singlets given by Eq.(4c). 

The ■P-projected dopons describe the physical doped 
carriers. Calculating explicitly the 'P-projected dopon 
operators we get 

Hu = Y.T ^^^l " ^S'^f) (i - 2S,r) d, 



(19) 



The representation (fTO|) has been used by Ribeiro and 
Wen within the mean-field approximation. However, 



they oversimplified the magnetic term in the following 
way, 

jJ2iSr ■ S, - i), (20) 



where J = J(l — x)^, x is a density of dopons. It is clear 
that the representation (j20p totally ignores the dynami- 
cally induced doping changes in the underlying spin cor- 
relations. The authors instead introduce phenomcnologi- 
cal x-dependent hopping parameters to take into account 
the feedback of the dopons on the spin dymanics. To take 
into account the actual J dependent spin-dopon interac- 
tion, one should use the constraint Ti = 0. It is satisfied 
provided Mi + fiiSi = 0. This yields (up to unessential 
constant factors). 



Hi 



tj 



2 '^2 

ija 



1 



2Sjf) dj 



[0^ ■ Sj - 1/4) + {S,Mj + SjA%) 



-{M,Mj - fiifij/A) 



(21) 



The magnetic term in this representation explicitly ac- 
counts for the spin-dopon interaction as produced by the 
magnetic moment-moment interactions. In general, it is 
the representation ((2T|) that should be used as a starting 
point to apply a mean-field approximation. In this way 
a complete dynamical mean-field phase diagram to de- 
scribe hole/electron doped cuprates emerges in contrast 
to a semi-phenomenological one previously derived^' 

Close to half-filling an alternative way of dealing with 
the constraint can be proposed. First, we can drop the 
projection operator and add the constraint with an ap- 
propriate Lagrange multiplier. That is, we now write: 



Hi 



tj 



d' d 



J 5^(5. • S, nf)il - nf) + At, (22) 



where A is to be send to -|-oo at the end of calculations. 
Second, since the dynamics is now restricted to the phys- 
ical subspace, we can close to half-filling make the change 
J — > J to get 

Hu J2 '^t'A.d,. + jY.{S,-S,-\) + XT. (23) 

ij<7 ij 

By writing the Hamiltonian in this way, we see the de- 
cisive role played by the constraint which incorporates 
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the interaction between dopons and lattice spins, while 
the non-constrained Hamiltonian simply corresponds to 
non-interacting dopons and lattice spins. 

The large-A limit eliminates unphysical states with the 
total spin j — I and at the same time dynamically gen- 
erates a vacancy on a lattice site in the form of a spin- 
dopon singlet-^. An analogy of this result can be drawn 
with the dynamical formation of the Zhang-Rice singlet 
produced by the hybridization effects that strongly bind 
a hole and a Cu^"*" ion together to form a local singlet 
state associated with such a vacancy^. The important 
point here being again the fact that such an empty site 
(vacancy) can be physically interpreted as a spin singlet. 

We now introduce the chemical potential for dopons 
and use the explicit representation of the constraint, Eq. 
(O, to get the Hamiltonian: 

ijcr ij i 

(24) 

where tl^^-^ = 2tij + (|A — ^i)5ij. The parameter A must 
be sent to +00 at the end of calculations^!. 

Finally, we can safely treat the constraint close to half- 
filling at the mean-field level. In this case, A is deter- 
mined sclf-consistently from the ground-state average 

i 

and becomes doping dependent. The mean-field Hamil- 
tonian obtained in this way is that of a Kondo-Heisenbcrg 
lattice problem, where the lattice spins play the role of 
localized magnetic moments while the dopons take the 
role of conduction electrons. 

In order to take into account a possible hybridization 
between the localized spins and the dopons^ one should 
use the following form of the constraint 

i i 

The relevant order parameter takes the form a :=< D >, 
where we have linearized the on-site operator product in 
the averaged constraint in the following way, DjOi w 

d] < D > + < > A- < >< D > . This 
procedure preserves the SU(2) symmetry, which must be 
present in the underdoped phase. Since a cx ^/x, the 
error produced by this simplification is at most of order 
0{x),x —> 0, which does not affect the results quanti- 
tatively in that regime. The breakdown of the Kondo 
regime implies then ct = 0. 

IV. PHYSICAL IMPLICATIONS 

In the previous section we showed that the t-J model 
can be mapped onto a Kondo-Heisenberg model for do- 
pons and lattice spins. In this section we explore some 



immediate physical implications of that mapping, leaving 
a more detailed analysis to a future work. 

The Kondo-Heisenberg lattice model has attracted 
much attention in the context of heavy-fermion 
systems^^i^liS^i^. It is believed that a variety of phys- 
ical phenomena could be accounted for by that model, 
such as non-Fermi liquid behavior, different types of 
magnetic, charge ordering and perhaps unconventional 
superconductivitj^. 

In fact, recent experiments have revealed strik- 
ing similarities between quasi two-dimensional heavy 
fermion systems (the CeMIns family) and the high- Tc 
cuprates^i^i^. The mapping discussed in the last sec- 
tion suggests that these similarities could be accounted 
for by the fact that both quasi-2(i heavy fermions and 
high-Tc cuprates capture universal features of strongly 
correlated electron systems in the presence of strong AF 
correlations. 

Our mapping reinforces earlier suggestions of a com- 
mon magnetically mediated mechanism for supercon- 
ductivity in heavy fermion compounds and in the 
cuprates24'-^^. In this way, the superconductivity in the 
cuprates can be directly associated with the pairing of 
dopons induced by the Kondo like interaction with the 
lattice spins. However, if this is indeed the case, a cru- 
cial question arises naturally: why are the critical su- 
perconducting transition temperatures observed in heavy 
fermions (low Tc) and cuprates (high Tc) so different 
from each other? Our explanation for that is signalled 
by the different charge carrier mass renormalizations and 
the typical coupling constant magnitudes in those two 
referred systems. Suppose the critical superconducting 
temperature is given generically by Tc = A exp (—cm* / A) 
where c is a constant, m* is the charge carrier effec- 
tive mass, A is the lattice spin-dopon coupling and A is 
some typical energy scale. In heavy fermion compounds 
A = Jk is small, and m* is two or three orders of magni- 
tude bigger than the bare electron mass. This leads to a 
very small Tc. In contrast, from infrared Hall measure- 
ments on underdoped LSCO and YBCO, m* is of the 
same order of magnitude as bare electron mass and A is 
large. This is due to the fact that in the large A limit, 
for optimally doped cuprates, we arrive at a Kondo like 
regime with Tc oc Tk oc D, the dopon bandwidth. 

Let us now connect our result more directly to the re- 
cent experimental results of Doiron-Leyraud et al^. Right 
at half-filling and below Neel temperature T^r oc J/4, the 
antiferromagnetic ordering is accounted by the Heisen- 
berg interaction term in Eq. ([24|) . Above T^r, thermal 
fluctuations destroy the long-range order. However, since 
the spin exchange energy J is in fact extraordinarily large 
the system still shows strong short-range AF correlations 
well above Tn- This phase is accounted for by the spin 
liquid state of the spin-spin singlets. As the dopon con- 
centration increases, the long-range AF order is melted 
by the quantum mechanical jiggling of the local spin mo- 
ments induced by the small (in this regime) dopon-spin 
interaction, and it eventually disappears altogether. Al- 
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though the RKKY spin-spin interaction induced by do- 
pons produces by itself the long-range AF ordering of the 
lattice spins, its strength is oc which is much less than 
the spin-spin exchange energy, J. As a result, at some fi- 
nite dopon concentration, the AF long-range order gives 
way to short-range AFM spin-spin correlations and the 
lattice spins become paired. As the doping increases, the 
individual lattice spins become less correlated with each 
other due to the competition between AFM fluctuations 
and the Kondo screening. 

Suppose we are now in the pseudogap regime. The lat- 
tice spins form singlet pairs interacting with the dopons 
by means of a still weak A coupling. The low density 
dopons are the only fermionic carriers present in the sys- 
tem which can be associated with the small hole pocket 
Fermi surface (FS) of the pseudogap state. The small 
volume of such a FS is accounted by the low density do- 
pons present in the system. As the density of dopons 
increases, the dopon-spin coupling also increases and the 
dopon-spin singlets evolve continuously out of the pseu- 
dogap state into a more Kondo like regime'^^. 

Let us now estimate the critical density associated with 
such a crossover. The necessary energy to break the lat- 
tice spin pairs is roughly J. Since the individual spins 
in each pair and dopons become closely coupled to each 
other by means of the increase of A , the dopon kinetic 
energy 2te, where x is dopon concentration, soon be- 
comes of the order of J. Consequently, when 2tx > J, 
the spin gap is destroyed and the resulting FS is now en- 
larged by the presence of the highly correlated charged 
spins which together with the dopons are now counted 
as charge carriers. Taking simply J/t « 1/3, we arrive 
immediately at the lower bound estimate for the criti- 
cal density Xc ~ 0.17, which is in very good agreement 
with the experimental value for the small versus large FS 
crossover which, for the hole doped cuprates takes place 
in the doping range of 0.15 — 0.25, and it is associated 
with the complete disappearance of the pseudogap state. 

It is also worthwhile at this stage to compare our 
scheme with the earlier mean-field slave-boson formu- 
lation. In that representation, the electron operator is 
decoupled into a spinon (fermion) and a holon (boson). 
Clearly when the spinous are paired into singlet states 
and the charged bosons are not yet Bose-Einstein con- 
densed the resulting spin gapped state has no FS to 
be associated with. This is in direct disagreement with 
the recent experimental findings which demonstrate the 
metallic character of the pseudogap phase. 

To complete the overall discussion of the physical im- 
plications of our dopon-lattice spin system, we need to 
clarify the onset and disappearance of the superconduct- 
ing state in both underdoped and overdoped regimes. 
In the underdoped limit, as emphasized earlier, the lat- 
tice spins form singlet pairs interacting weakly with the 
dopons through A. This interaction A naturally leads 
dopons to condense into BCS like pairs at temperature 
T < Tc- Notice that since the coupling A is weak at 
low doping, the superconducting gap resulting from BCS 



condensation of dopons is strongly doping dependent and 
it is not directly related to the spin gap. As a result the 
superconducting gap and the spin gap (pseudogap) are, 
in practice, two independent energy scales at very low 
doping. This is in agreement with recent angle-resolved 
photoemission (ARPES) and Raman experiments which 
distinguish the roles of the nodal and anti-nodal gaps 
in the low doping superconducting phase. In contrast, in 
the optimal doping region for larger dopon concentration, 
the NDO constraint must be treated with care and as a 
result A grows accordingly. Such a growth of A strongly 
ties the spin and dopon to each other. Therefore, at suf- 
ficiently large dopings, the spin and dopon gaps should 
become indistinguishable from each other. The supercon- 
ducting phase is well described in this regime by a single 
energy scale. As a result at a sufficiently large A, both 
gaps are destroyed simultaneously and we end up with a 
low-energy Kondo like metal^I-. Needless to say, the pro- 
jection NDO constraint is a crucial ingredient in all our 
arguments and it allows us to give a simple explanation 
of important recent results. 



V. CONCLUSION 

In the present paper we discuss the physical meaning 
and some implications of the theory of the projection 
constraint in the doped carrier representation of the t-J 
model. The basic conclusions that can be drawn from our 
consideration are as follows. Firstly, the complete theory 
that incorporates the constraint sounds quite appealing, 
since it allows for a controlled mean-field treatment of 
the t-J model in the most interesting region close to half- 
filling. This happens at the expense of having more com- 
plicated interaction terms in the spin-dopon t-J Hamil- 
tonian. It would be of utmost importance to derive the 
complete mean-field theory in this representation which 
is now under consideration. However, as a first step in 
that direction we just discuss some immediate qualitative 
consequences of the improved spin-dopon approach. 

Secondly, the constraint enforced by the Lagrange mul- 
tiplier term allows for an explicit mapping of the t-J 
model onto the Kondo-Heisenberg lattice model in the 
underdoped region . This indicates that the physics be- 
hind these two model are indeed related to each other. 
This mapping is very appealing in view of recent experi- 
ments that suggest striking similarities between quasi two 
dimensional heavy fermion systems and high- Tc cuprates. 
Some physical implications are briefly discussed, point- 
ing to the unified physics of heavy fermion and high- 
Tc materials. Namely, it is possible that the very same 
physical mechanism is responsible for the formation of 
Cooper pairs in those systems, with different critical tem- 
peratures related to different mass renormalization of the 
charge carriers and to different magnitudes of the existing 
coupling constants. 

We also discuss the small-large evolution of the FS with 
doping. We associate this crossover to the closure of the 
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spin gap and the destruction of the pseudogap state. We 
estimate the lower bound density for such a crossover. 
We make direct contact with the recent FS experiment 
of Doiron-Leyraud et al. When the spin gap is present, 
the pseudogap state has a small FS with the dopons be- 
ing the only available charge carriers. With the increase 
of the dopon concentration and with the consequent in- 
crease of the coupling A, the dopons and the individual 
lattice spins become strongly correlated to each other. In 
this way when the spin gap is destroyed, above the su- 
perconducting temperature T^, both dopons and lattice 
spins generate the associated FS. 

Making explicit use of the projection constraint we also 
discuss the onset and disappearance of superconductiv- 
ity in both underdoped and overdoped regimes. Further 
work is however needed to explore in a more quantitative 
basis our new mean-field doped carrier formulation of the 
t-J model. This is already in progress and the results will 
be presented elsewhere. 
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VII. APPENDIX: ANALOGY WITH THE 
GUTZWILLER PROJECTION 

The global projection operator described above is sim- 
ilar to that used in the U = oo Hubbard model. Let 
us consider the Hamiltonian Huub = J^ija^ij'^lcr'^j'^ + 
^X]i"«T"«i- ^ ~* system is sub- 

ject to the constraint -I- rii i < 1. This constraint is 

equivalent to Tf = rn^nii — 0. In this way when Tf 
acts on the unphysical state (doubly occupied) we have 
Tf\unphys)i = \unphys)i. Therefore, = 1 — Ui^riii 
is a projection operator that eliminates the unphysical 
state at site i. The gauge transformation generated by 
this constraint, 

leaves the projected electron operators = Ca{l — 
ri-a-) intact. The global projection operator is the 
well known Gutzwiller projector P*^ = IliP[^. We 
can then implement the constraint writing Hnub ~ 
P^J2tjcr^v4cr(^j<yP'^' which is equivalent to Hnub = 

X^ijtT ^ij^ia^jo- From this point of view, the parameter U 
of the Hubbard model becomes the Lagrange multiplier 
which is necessary to enforce the NDO constraint. 
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